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Abstract 

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating 
Boussinesq equations are analyzed in the asymptotic limit of strong stratification in the Burger number 
of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It 
is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the 
full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with ver- 
tical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes 
on inertial-gravity waves are evidenced. The ‘split’ of the energy transfer of the vortical and the wave 
components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero 
to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy 
spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the 
anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. In- 
termediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where 
the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) 
of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical 
viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and 
geostrophic departure) are obtained. 
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1 Introduction 


The turbulent flows that are subject to rotation and stratification have many important applications in 
geophysics and engineering (Fernando and Hunt, 1996; Hopfinger, 1989). An important class of geophysical 
flows can be characterized as strongly rotating and strongly stratified with both effects playing an important 
role in the dynamics. This is the so called Burger number of order one regimes where the effects of rotation 
and stratification enter at the same order in asymptotics (McWilliams 1985). One of the major difficulties 
encountered in understanding dynamics of geophysical flows is the influence of the oscillations generated by 
the rotation and stratification. One major effect of rotation and stratification is through “phase scrambling” 
for the wave phase. In this paper rotating Boussinesq equations are analyzed in the asymptotic limit of 
strong stratification in the Burger number of order one situation as well as in the asymptotic regime of 
strong stratification and weak rotation. 

The very useful and thought provoking multi-scale analyses of rotating/stratified turbulence is presented 
in Riley et al. (1981), Lilly (1983), McWilliams (1985). In particular, they argue that the velocity field 
of a rotating, stably stratified fluid may be regarded as a superposition of waves which are modulated on 
the longer turbulence time scale. In our approach, the collective contribution to the dynamics made by 
waves is accounted for by rigorous KAM (Kolmogorov-Arnold-Moser) type theory and rigorous estimates 
of wave resonances and quasi-resonances via small divisors analysis. Our theory handles rigorously 3-wa.ves 
resonances, but goes much deeper into the structure of quasi 3-wave resonances and their contributions. 
This mathematical approach was initiated in Babin, Mahalov and Nicolaenko (henceforth BMN) (1995), 
Mahalov and Marcus (1995) in the context of geophysical flows. In Bartello (1995), the relative physical 
importance of different resonances is discussed in depth. In this paper we present the physical predictions 
and implications of our rigorous mathematical analysis. Interactions between internal waves and the vortical 
(quasi-geostrophic) modes remain as one of the important questions to be addressed by strong interaction 
theory (Muller, Holloway et al ., 1986; Warn, 1986; Farge and Sadourny, 1989; Lelong and Riley, 1991). 

The governing flow equations for rotating stably stratified fluids under the Boussinesq approximation are 

<9 t U -f U • VU + 2f2 0 e3 xU = -Vp+/>i e 3 , V • U = 0, (1.1) 

dtPi + V-V Pl = -NgU 3 , (1.2) 

where rotation and mean stratification gradient are aligned parallel to e 3 = [0, 0, 1]. Here U = (U 1 , U 2 , U 3 ) is 
the velocity field and p\ is the buoyancy variable; No is the Brunt- Vaisala frequency for constant stratification 
and il 0 is the frequency of background rotation, /o = 2fio- We focus on inviscid Eqs. (1.1)-(1.2) or with 
small uniform viscosities. 
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Now we introduce useful non-dimensional parameters. Let Uh be a characteristic horizontal velocity 
scale. Let H and L be vertical and horizontal length scales and a 3 = H / L is the aspect ratio parameter. We 
define Froude numbers based on horizontal and vertical scales: 

F h = U h /LN 0 = l/N, F v = U h /HNo = F h /a 3 . (1.3) 

The classical Rossby and anisotropic Rossby number are defined as follows 

Ro = U h /2Lti o = 1/20 = 1//, Ro a = a 3 Ro. (1.4) 

The time is dimensionalized using the turbulence time scale L/Uh • In the Burger 0(1) regime Ro a ~ F h . 

We are not taking a 3 — > 0; rather its value is fixed by shallowness of the atmosphere, a 3 ~ 5 x 10 -3 to 10 -2 
for synoptic scales. Its smallness effectively downsizes Ro a ~ 5 x 10~ 4 to 10~ 3 , as Ro ~ 0.1 in midlatitudes. 

For mesoscales, L ~ lOOfcm, a 3 ~ 10” 1 and Ro a ~ 10 -2 . The anisotropic Rossby number Ro a or/and the 
vertical Froude number F v control our uniform error estimates. 

The Burger number characterizes relative importance of the effects of rotation and stratification (McWilliams, 
1985): 

Bu = Rol/FS = Ro 2 / F 2 = N 2 a 2 /f 2 = N 2 a 2 /f 2 . (1.5) 

In Eqs. (1.3)-(1.4) / = Ro -i and N = Ffl 1 are dimensionless rotation and Brunt-Vaisala. parameters, 
respectively. The relative importance of rotation/stratification is measured by the Burger number with 
Bu « 1 corresponding to rotation dominated and Bu » 1 corresponding to stratification dominated 
flows. Herring and Metais (1989) observed horizontal layering of the velocity field in numerical simulations 
of stratified turbulence, while Bartello et al. (1994) noted the formation of quasi- two-dimensional structures 
in rotating turbulence. 

Fourier series will be used in this paper to represent physical fields in a parallepiped [0,27r] x [0,27r/ao] x 
[0,27r/a3], 0 < a -2 < 1. We denote k^/a.n by k 2 , k 3 /a. 3 by k 3 . Following Bartello (1995), it is useful 
to distinguish between three sets of wavevectors k = (&i, k 2 /a 2) k 3 /a 3 ), |Ar| 2 = k 2 -f k\/a\ + k\ja\\ the 
barotropic set {k : k 3 = 0), the set with only vertical variability {k : ki = k 2 — 0}, and the remaining 
baroclinic vectors {k : k 2 + ki ^ 0 , A ?3 ^ 0). Then the operation of vertical averaging corresponds to 
projection on barotropic fields; the operation of horizontal averaging corresponds to projection on fields with 
only vertical variability. In this paper as well as in our previous work (BMN, 1996a, b, cl) we emphasize the 
important role of operations of vertical and horizontal averaging in investigations of rotating and stratified 
flows. 

Regimes of geophysical dynamics presenting the global picture for small Froude or Rossby numbers are 
shown in Figure 1 which summarizes the physical implications of our mathematical analysis. Since we are 
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not taking a .3 — >• 0, either Fh or F v can be used in description of asymptotic regimes. Then Fr denotes 
either of these numbers. The rotation dominated case corresponding to Ro a — » 0, Bu — > 0 and Fr finite 
was considered in BMN (1995), Zhou (1995), BMN (1996a, c, d, e), Mahalov and Zhou (1996) (Figure 
1, vertical axis). In this case we proved the generalized Taylor-Proudman theorem establishing splitting 
between vertical averages of U, p\ (two dimensional- four component, 2D-4C barotropic fields) and reduced 
ageostrophic field. Following Reynolds and Kassinos (1996) 2D-4C refers to fields with four components 
depending on two variables aq and a? 2 ; 3D-2C refers to fields with two components depending on three 
variables xi, x 2 and X 3 etc. By splitting we mean that the barotropic field decouples from the ageostrophic 
one, without feedback from the latter onto the former, to the lowest order. In the limit AT — > 0, the usual 
quasi-geostrophic field does reduce to vertically averaged fields; this is what is meant by “geostrophic” in this 
limiting context. In this limit vertical shearing is fully controlled which is reflected in adiabatic invariants 
associated with vertical shearing (exact conservation laws in reduced equations), with = (U,pi) 


dt 


dU f 

dx?, 


}j dxidxodx3 = 0 . 


( 1 . 6 ) 


In the limit N — > 0 the temperature decouples from the dynamics and behaves as a passive scalar. The 
dynamics of vertically averaged velocity fields reduce to classical 2D-3C Euler equations and are subject 
to inverse energy cascades in U 1 , U 2 as in 2D turbulence. For Bit — > 0 we have shown rigorously that 
energy cascades for the ageostrophic (AG) field are completely frozen in X 3 and the dynamics is pure phase 
turbulence (BMN, 1996a, c,d,e); freezing of energy cascades in X 3 for the “baroclinic” component follows 
from (1.6). In pure phase turbulence, the amplitudes of the ageostrophic modes remain constant in absolute 
values; turbulent dynamics are restricted to the phases of the ageostrophic modes. The ageostrophic field is 
phase locked to phases associated with vertically averaged vertical vorticit.v and vertical velocity which are 
advected by 2D turbulence of vertically averaged fields. There is no slaving of the amplitudes of ageostrophic 
modes by the 2D turbulence, only phase locking. We calculated Doppler phase shifts induced by turbulence 
of vortical (vertically averaged) fields on inertio-gravity waves in this limit. In the case of 3D rotating Euler in 
the small anisotropic Rossby number situation we described regimes with no energy flux in the ageostrophic 
(AG) component and formation of K AM-type regimes with frozen in X 3 ageostrophic cascades (pure phase 
turbulence, frozen turbulence). Similar freezing of energy cascades was observed by Farge (1988) in the 
context of rotating shallow- water equations and in Pushkarev and Zakharov (1996) in numerical experiments 
describing turbulence of capillary waves. 

Next is the regime of strong rotation and weak stratification as shown in Figure 1. Besides the operation 
of vertical averaging there is a piece of 3DQG (quasi-geostrophic, Pedlosky 1987) which plays an important 
role in the dynamics. It is formally obtained by expanding 3DQG equations in a small parameter N/f. 
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This procedure is similar to the one described in Section 4 for the regime of strong stratification and weak 
rotation. The corresponding reduced equations, higher order corrections and mathematical convergence 
results for the case of balanced and unbalanced initial data are presented in Avrin, Babin, Mahalov and 
Nicolaenko (henceforth ABMN) (1996). 

As the effects of stratification are increased further (see Fig. 1) AG cascades become possible. In the 
limit of strong rotation and strong stratification corresponding to Ro a — >• 0 , Fr — > 0, but Bu = 0(1) 
we established splitting between 3DQG and the reduced ageostrophic field (BMN, 1996b) confirming the 
similar splitting for rotating shallow- water equations which we obtained in Mahalov and Marcus (1995) 
and for rotation-dominated case in BMN (1995), (1996a). Again, by splitting we mean that the QG field 
decouples from the ageostrophic one. Djmamics of AG further splits along uncoupled resonant quadruplets of 
Fourier rays with AG energy conserved on each resonant quadruplet of rays. Energy cascades are now allowed 
(unfrozen) for the ageostrophic field but they are restricted to families of 4 rays in Fourier space. Energy 
transfers (direct and inverse cascades) are not prohibited but restricted on uncoupled families of resonant 
four rays for AG field. Direct energy cascades of AG field provide mechanism for nonlinear geostrophic 
adjustment. This is fundamentally different from the rotation dominated regimes where AG cascades are 
frozen. Spectral differential molecular viscosities for QG and AG fields with explicit dependence on the 
rotation/stratification parameters behave differently. Let zq and zq be the kinematic viscosity and the heat 
conductivity, respectively; the ratio zq/zq is known as the Prandtl number. Through a simple computation 
of 2-wave resonances in the Craya-Herring basis (Lesieur, 1987) the effective differential spectral molecular 
viscosities vqg and vag are given by 


m2 


"««(") = 1'2 + (H - Vl) ||V p + y . •'AG(n) = ^ + (v 2 - ^) — 


(1.7) 

/rug ' |n'p-fp-ng 

where p = f /N, |tz/| 2 = ?zj n\/a\. It shows that dissipation affects QG and AG fields differently. This 
impacts on direct numerical simulations of QG fields in the context of numerical simulations of atmospheric 


flows. 


Partial control of vertical shearing is obtained allowing us to prove long time existence of solutions of 
inviscid Boussinesq equations (BMN, 1996b and Section 6). Also, a flow which is initially wave dominated 
remains wave dominated even through decay (confirming Metais et. al, 1996). In Section 3 we show that 
horizontally averaged buoyancy variable pi is an adiabatic invariant (this result is true for all resonances 
including 3 wave resonances); we calculate Doppler phase corrections associated with to a linear profile 

(constant iV 0 ). This confirms and generalizes the work of Legras (1980) on phase shifts who showed the 
existence of statistical dephasing effects induced by turbulent processes on Rossby waves. Frequency shifts 
induced by turbulent processes on inertial waves were calculated in BMN (1996a); the case of frequency 
shifts induced on waves in rotating shallow- water equations was considered in BMN (1996c, e). As the 
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effects of stratification are increased (see Fig. 1) vertical shearing dynamics in the ageostrophic field increase 
and is conveniently characterized using the divergent velocity potential x(U ®i, X 2 , £ 3 ) given by the formula 
X — {—Ah)~ 1 dU 3 /dx 3 which is coupled to geostrophic departure. The geostrophic departure —V\p\ + 
f-£rCurW • e 3 characterizes imbalance in the vertical motion or omega equation (Eqs. ( 1 . 16)- (1.17)). Up 
to a normalization, the divergent velocity potential y, the geostrophic departure and 3DQG mode form the 
Craya-Herring cyclic basis which is used in this paper to represent physical fields. This is further described 
below. 

In the inviscid regime Fr — > 0, Ro a = 0(1), Bu — > +00 (Figure 1, horizontal axis) we prove that there is 
no bound on vertical shearing associated with the dynamics of 3D-2C (3 dimensional, 2 components) decou- 
pled pancakes (parametrized in £ 3 ) with different pressures at every level; this leads to unbalanced dynamics 
at the lowest order. There is no saturation of the exponential build-up of vertical enstrophy (in small vertical 
scales) for AG dynamics as the latter is coupled to the quasi 2D field thru dUq 2D /dx 3, dUq 2D fdx 3. The 
major problem is lack of boundedness of vertical shearing in quasi 2D equations (Lilly, 1983). We show that 
horizontally averaged XJ 1 (x 3 ), U 2 (x 3 ), pi(x 3 ) are adiabatic invariants providing a feedback onto AG turbu- 
lence. U 1 , U 2 are adiabatic invariants only if Q = 0; otherwise, they undergo rigid 0- rotation (this result 
holds for all resonances including the 3 wave resonances). However, these adiabatic invariants are not enough 
to saturate vertical shearing. Worse the lack of boundedness of dUq 2D /dx 3 , dUq 2D /dx 3 leads to explosive 
exponential growth of the AG dynamics. Of course, control of vertical shearing can be achieved trivially by 
introducing vertical viscosity; however, this corresponds to a non-physical laboratory set-up rather than the 
real atmosphere (A. Majda, private communication), or a poorly-resolved (in £3- scales) numerical model 
(P. Bartello, private communication). 

In the intermediate asymptotic regime corresponding to strong stratification and weak rotation (Bu — * 
+ 00 , f /N small) the effects of weak rotation on the dynamics are accounted for by an asymptotic expansion 
in a small parameter p — f /N (Section 4). Full saturation of vertical shearing is obtained for all times. 
Equations describing balanced dynamics are intermediate between 3DQG equations valid in the regime Bu = 
0(1) and qua.si-2D decoupled pancakes without any control of vertical shearing (Lilly, 1983). In this paper we 
show how weak rotation regularizes vertical shearing and calculate its effects on AG dynamics. Our reduced 
equations have a conservation law associated with vertical shearing which allows to control AG vertical 
scales for all times. There is no need to resort to vertical viscosity as the principal stabilizing mechanism 
(Reynolds number Re ~ 10 12 in atmospheric flows). Two scalar decoupled equations for ageostrophic 
components (divergence velocity potential and geostrophic departure) are obtained. These equations have 
coefficients depending on regularized quasi-2D fields and can be used for more accurate and robust numerical 
simulations of geophysical flows in the regime of strong stratification and weak rotation. AG dynamics is 
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driven by regularized vertical shearing of the pancakes. 

In this paper we emphasize physical predictions which follow from rigorous mathematical analysis of 
Boussinesq equations in the strongly rotating/stratified Bu = 0(1) regime as well as in the asymptotic 
regime of strong stratification and weak rotation. The mathematical theory is based on rigorous small divisor 
estimates and KAM type (Kolmogorov-Arnold-Moser) theoretical considerations to rigorously control wave 
resonances, especially the 3- waves of the fast-fast-fast resonances (Figure 2). 

On the physical side, for the Bu = 0(1) regimes discussed in Sections 3 and 5, we establish statistical 
dephasing effects induced by turbulent processes on inertio-gravity waves with 3DQG turbulence acting to 
renormalize both frequency and viscosity of the waves. We generalize the work of Legras (1980), Carnevale 
and Martin (1982). In particular, we calculate Doppler phase corrections associated with the fact that 
horizontally averaged buoyancy variable p\ is an adiabatic invariant. Namely, pi{xz) has a O(e) variation on 
large times when Ro a ~ Fh ~ e. Rigorous mathematical analysis based on small divisor estimates shows that 
3 waves fast-fast-fast resonances are rare in the Bu = 0(1) atmospheric-jegimes (BMN, 1996b and Section 
6) as well as in Bu « 1 regimes (BMN, 1996a). In fact, just switching on even weak rotation destroys 
the 3- waves resonances found in the pure stratified case / = 0. One of the hardest parts of our analysis is 
to estimate the total probabilities of quasi-resonances, that is the width of Arnold tongues coming out of 
points (set of measure zero) where 3 wave resonances are possible (see Figure 2). These resonances are not 
neglected but rather weights are assigned to them according to their importance (BMN, 1996d and Section 
6). Even 3 waves resonances do not alter the global picture: they correspond to higher order corrections 
(“Arnold drift”). 

The width of Arnold tongues scales algebraically with Ro a (see Corollary 6.1 and Figure 2). For synoptic 
scales at midlatitudes, L « 2000 km, H ss 10 km, the width (normalized probability of a 3 wave resonance) 
is of order (H/L) 2 = (5 x 10~ 3 ) 2 and it is of order (10 -2 ) 2 at mesoscales ( L « 100 km); 3 wave interactions 
become significant only for tall columns (see Section 6). In regions free of fast-fast-fast interactions (set of 
full Lebesgue measure) nonlinear geostrophic adjustment takes place via “catalytic” interactions between one 
QG mode and two AG modes confirming the insight and numerical simulations of Bartello (1995). Bartello 
(1995) also discusses fast-slow-slow interactions, which are non-resonant. In our work, they appear at the 
next order in Ro a or 1/N at Bu = 0(1), and contribute to the feedback of the ageost.rophic field onto 
the QG one (ABMN, 1996). Inside Arnold tongues where fast-fast-fast 3 wave interactions are possible we 
expect Arnold drifts associated with neglected higher order resonances. Guided by KAM theory in finite 
dimensional systems and the fact that Arnold tongues are very narrow in our problem (Section 6), we expect 
that these drifts will be slowly evolving (cf. Figure 2). Thus we expect the dynamical picture to be intact even 
inside Arnold tongues where 3 wave interactions are possible. Our analysis (Section 6) proves that 3 wave 
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interactions will not contaminate large scale dynamics. Breakdowns can occur only locally corresponding to 
small horizontal scales. In this case a local anisotropic Rossby number based on local horizontal scale will 
become large and both local breakdown and collapse leading to fully 3D dynamics (locally) will be possible. 
We have analyzed such 3D instabilities in an idealized case of rotating columnar flows (Mahalov, 1993). Here 
our rigorous mathematical analysis amounts to the following: the probability of such a localized breakdown 
to 3D turbulence to extend to larger scales via a catastrophic inverse cascade is very small. This is indeed 
what is observed in the atmosphere, where near balanced dynamics are sustained for all times in the larger 
scales. The regime Fr — » 0, Ro a = 0(1), Bu — > oo requires a special attention. In this regime 3DQG 
equations degenerate to quasi-2D equations lacking control of vertical shearing and leading to unbalanced 
dynamics at the lowest order. To this end, a regularization based on weak rotation is proposed in Section 4. 

In BMN (1996b) we established splitting between 3DQG and reduced ageostrophic field; here we show 
that the structure of reduced ageostrophic equations (via. Craya-Herring cyclic basis) implies upscale (inverse) 
energy transfer of rotational (vortical) energy via. QG mode versus downscale cascade of wave energy via 
slow-fast-fast. resonant catalytic interactions (following Bartello, 1995 notations). In cases Bu = 0(1) and 
Bu >> 1, AG field satisfies uncoupled families of equations on 4 rays in Fourier space. As the effects of 
stratification are increased direct cascade of AG energy (along these rays) toward small scales is enhanced. 
In Section 5 we introduce the anisotropic phase coherence tensor and model anisotropy in Bu — 0(1) regimes 
of geophysical turbulence. There is a spectral gap (i.e. a power law scaling break) between QG and AG 
spectrum with AG spectrum being shallower than the typical k~ 3 QG spectrum; it varies smoothly between 
k~ 2 and ft -5 / 3 which is in agreement with numerical simulations of Metais et al. (1996), Bartello et al. 
(1996). 

A possible emerging picture of Burger 0(1) turbulence is that 3DQG turbulence being a guiding center is 
corrupted by phase turbulence and Doppler phase shifts; with the dynamics of the AG field being constrained 
to a uncoupled 4-rays families, with direct AG cascade restricted to the latter. The feedback of AG field on 
QG dynamics can be computed at next order in Fr or Ro a (ABMN, 1996). 

On the mathematical side, the challenge is to prove that dynamics of limit equations uniformly closely 
approximate geophysical turbulence in regimes: 

1. N = 0(1), Q >> 1 (rotation dominated regime), Fr = 1/N . 

2. A r >> 1, 0 >> 1, Bu = 0(1) (strongly rotating/stratified Burger number of order one regime). 

One needs uniform error estimates in Froude and Rossby numbers. It is important to realize that “weak” 
convergence results, “filtered” and o(l) error estimates cannot achieve this. In BMN (1996 a,b,d) and here 
(Section 6) we obtain strong convergence results with uniform error estimates on sets of full Lebesgue measure 
with initial data being in the Sobolev space H‘ (and even H 4 ) for both Bu = 0(1) and rotation dominated 
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regimes. This is in contrast with the work of Embid and Majda (1996) where only local existence results are 
established. Following general theorems of Schochet (1994), they state a pointwise convergence theorem with 
o(l) error on a time interval [0,T]; since they include (for Bu — 0(1)) 3-wave resonances they cannot give 
any uniform error estimates in the parameters. Our theory handles rigorously 3-waves resonances, but goes 
much deeper into the structure of quasi 3-wave resonances and their contributions. Moreover, the interval 
[0,T] in Embid and Majda (1996) is small where T can be as bad as any local time of existence of original 
3D Euler/ Boussinesq equations , and it is in fact limited by potential 3D vortex singularities of the full 3D 
Euler. Embid and Majda (1996) miss the regularizing effects of strong rotation and stratification. As Ro a 
and Fh -A- 0, Bu — 0(1) we prove that T — » +oo for (1.1)-(1.2). We agree with them only as to the QG 
operator decoupling, which we obtained independently in BMN (1996b) for Boussinesq equations. 

Now we describe the structure of reduced equations which will be derived in Section 3 for Bu = 0(1). 
In the reduced equations the total velocity splits into the quasi-geostrophic field wqg(^) satisfying 3DQG 
(quasi-geostrophic) equations 


^wq G = B 1 (wq G ',wq G ), (1.8) 

and the ageostrophic component satisfying in principle general equations of the type: 

dtWAG = B 2 (wq G ,W j 4 G ) -1- B 3 (W J 4G, W^g)* (1-9) 

In the case when the balanced QG dynamics are absent, w qg{1) — 0, Eqs. (1.9) reduce to equations describing 
wave turbulence (e.g. Zakharov et ai, 1992): 

<9 f .w AG = B 3 (w^ G ', w^ g ). (1-10) 

We prove that (1.9) holds only within very narrow Arnold tongues (Figure 2) with infinitesimal probability 
(Section 6); we do not just drop B 3 with some ad-hoc scaling argument. In the real atmosphere situation a 
near balance state is sustained and the quasi-geostrophic (balanced) field is present and it is dominant. In 
Babin, Mahalov and Nicolaenko (1996b) and here (Section 6) it is shown using small divisor estimates that 
for almost all Bu and a 3 = H / L only “catalytic” interactions rule AG dynamics: 

d t v?AG = B 2 (wqg, w^g) (1.11) 

where w qg{1) is a solution of 3DQG equations. Such interactions do not influence the slow QG modes, but 
act to transfer fast AG energy downscale (Bartello, 1995). Further, for all Bu, all a, 3 and almost all domain 
aspect ratios 02 , Bo(wqg, w^g) splits in Fourier space into uncoupled, restricted interaction operators on 4 
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rays families 


■ m l \ 


( ±1 0 0 ^ 


n, ) 

m -2 

= A 

0 

If 

1 — * 
0 


n 2 

{ m 3 ) 


^ 0 0 A1 j 


V ns ) 


In Eqs. (1.11) direct cascades of energy are allowed for w ag through B2 (wqg,wag)- Wave energy cascades 
toward smaller scales and is subject to strong dissipation in AG mode, via the anisotropic viscosity uag 
in (1.7). 

An important observation on the nature of interactions in Eqs. (1.11) and their impact on AG dynamics 
(AG cascades/frozen AG turbulence) is now in order. The nature of interactions in Eqs. (1-11) is funda- 
mentally different in rotation dominated and stratification dominated regimes. In the rotation dominated 
case (see Fig. 1) corresponding to Ro a — » 0, Bu — > 0, and Fr finite interactions in Eqs. (1.11) for n- 
th, n — (ni,n2,«3) Fourier mode are restricted to wavenumbers m — (mi, m2, m3) satisfying m 3 = ±713, 
|m/| = 1 77/ 1 , k + m = n (BMN, 1996a). Thus, modes m which interact with a given mode n lie on the same 
energy shell as n. Interactions in (1.11) are localized in Fourier space and do not extend to zero or infinity. 
This special nature of interactions in the rotation dominated regime implies freezing of energy cascades in 
X 3 in Eqs. (1.11). Similar freezing of energy cascades of AG field (frozen AG turbulence) was proven for 
rapidly rotating shallow-water equations in Mahalov and Marcus (1995), BMN (1996c). In particular, it was 
shown that in the case of rapidly rotating shallow-water equations interactions in ?V = (ni, n2) AG equation 
are restricted to wavenumbers m' = (??ii, mf) satisfying \m'\ — |?V|, k' + m' — n' . 

The nature of interactions in (1.11) drastically changes when stratification is increased (see Figure 1). As 
shown in this paper, in the Bu = 0(1) as well as Bu » 1 regimes interactions for the rc-th mode in (1.11) 
are restricted to four rays families given by (1.12) for almost all 02 - Thus, interacting modes lie on families 
of four rays connecting 0 and infinity in Fourier space. AG energy cascades are now possible along these 
rays from small wavenumbers (large scales) to large wavenumbers (small scales). 

Vertical shearing operations in B2 (wqgs w^g') ai ' e conveniently expressed in terms of divergent velocity 
potential y. We recall that % and 4 } are related to a horizontal velocity field U/> by the formula U/, = 
£3 x VftV’ + VhX- The operator Bo(wQG,w i 4G) takes an especially simple form in the Craya-Herring cyclic 
basis. In this basis the first component of wag is simply the divergent velocity potential \ = (— A/J" 1 
or in Fourier space 


. h 3 1 ! 

Xn — * I - I I ~ /\ w n 


(1.13) 


where h' — (n 1 ^ 2 / 02 ), and the components of the full field in the Craya-Herring basis are (w° , uff). 
Here w® is the component of the quasi-geostrophic mode and uf , iif are the ageostrophic components. The 
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important role of x in our AG operator B 2 ties with the more classical theory of balanced models (Gent and 
McWilliams, 1983a and 1983b) and the findings of Farge and Sadourny (1989), Farge (1988) on effects of 
inertio-gravity waves and rotation on two-dimensional turbulence. 

The second component ur of w a.g in the Craya-FIerring basis is related in a simple way to the vertical 
motion or omega equation. We recall that the vertical velocity does not appear in 3DQG equations and the 
whole purpose of manipulations leading to the omega equation is to obtain a diagnostic equation from which 
vertical velocity component can be calculated (e.g. Holton, 1992; Gent and McWilliams, 1983a). We briefly 
recall these calculations starting from Eqs. (1.1)- (1.2) written in the form (/o = 2Qo) 

d t U + f 0 e 3 xU-pi e 3 = -Vp - U • VU, V • U = 0, (1.14) 

d t .pi + NqU 3 — — U -V p\ . (1.15) 

First, we apply the operator f 0 -^curl to Eq. (1.14) and obtain a scalar equation by projecting the result on 
e 3 . Then we apply horizontal Laplacian to Eq. (1.15). Recalling that curl(e 3 x U) • e 3 — we obtain 

by subtraction 

d,(-Vlp 1 + /o^u> 3 ) - (A r„M + = G (116) 

where G denote nonlinear terms; u > 3 = curl U • e 3 is vertical component of the vorticity vector. The omega 
equation is obtained by dropping d t (— +/o ^^ 3 ) in Eqs. (1.16). Then we obtain a diagnostic equation 
which is used to determine vertical velocity U 3 : 

-(Af 0 2 Vl + /„ 3 £|)t /3 = G. (1.17) 

Up to a normalization constant, second mode in the Craya-Herring basis is precisely the geostrophic departure 
— (buoyancy) + /o^d_ (vertical vorticity). It characterizes unbalance in Eqs. (1.16)-(1.17). The geostrophic 
departure has the form |n , | 2 /?i t ra — fon 3 (n x Ujj • e 3 = No\h'\ 2 p n — fon 3 (h x U|() • e 3 in Fourier space. Since 
N 0 \n'\ 2 p n -f 0 n 3 (h x Uj,) ■e 3 = WnlnUn'j ■q 2n , the geostrophic departure up to a normalization is precisely 
the mode q 2 in the Craya-Herring cyclic basis (see Eqs. (2.9)). 

This paper is organized as follows. In Sec. 2 we recast the Boussinesq equations in the Craya-Herring 
cyclic basis. This intrinsic coordinate system is particularly convenient for investigations of adiabatic con- 
servation laws and global decoupling. Next in Sec. 3, we develop a methodology and procedure for studying 
rotating Boussinesq equations in the strongly stratified limit. Here, the dynamical decoupling between the 
vortical and inertial-gravity wave components of the total field is achieved. It is shown that horizontally 
averaged buoyancy variable is an adiabatic invariant in the asymptotic state and phase corrections to buoy- 
ancy time scale which are associated with vertical shear of the horizontally averaged buoyancy are obtained. 
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Ill Section 4 we analyze an intermediate asymptotic regime of strong stratification and weak rotation with 
effective saturation of vertical shearing of pancakes by weak rotation (no vertical viscosity needed). The 
energy spectrum and spectral eddy viscosity models are deduced in Section 5. Finally, in Section 6 we 
describe the uniform convergence results and the regularization of Euler-Boussinesq equations, as Fh — » 0, 
Ro a — > 0, Bib = 0(1). These results require much less differentiability than in BMN (1996a): now only four 
derivatives on initial data are required for the most general convergence results on arbitrary long times. 


2 Boussinesq equations in the Craya-Herring cyclic basis 


We introduce a change of variables pi = N p (Metais and Herring 1989) and combine velocity and buoyancy 
variable in one variable = (U,p) after which Eqs. (1.1)-(1.2) written in non-dimensional variables take 
more symmetric form: 


<9 f U f + U • VU f + /RU 1 = -V ] p - N SU f , V • U = 0 
where V^p — (Vp, 0) and time was non-dimensionalized using turbulence time scale L/Uh ■ Here 


R = 



S = 



, J = 


0 


( 2 . 1 ) 


( 2 . 2 ) 


and R n , S n will denote the action of R and S on n-th Fourier component, M = / R + ATS. In this section 
we write Boussinesq equations in the Craya-Herring cyclic basis and use this representation to perform time 
averaging. Linear equations describing inertio-gravity waves are 


d t U f + MU 1 = -V f p, V • U = 0. 

Dispersion relation for inertio-gravity waves which are solutions of Eqs. (2.3) has the form 




2 / J 
n 3f a Z 


ko: 2 (bu^ + 


J\2 


N- 


+ Bu 


-l 


I nr I nr [n |“ 1 77 1 " |np 1 71 

;/|2 _ , „2 / ,,2 


5) 


(2.3) 


(2.4) 


where |t7| j = n\ + n^/ai + 77.3/05, \h'\ 2 — n{ -f n^j a\. Here Ro a is the anisotropic Rossby number and Bu is 
the Burger number defined by (1.4), (1.5). It follows from (2.4) that the effect of rotation and stratification 
are not uniform on scales. In the case | ?? / 1 / 1 713 1 >> 1 gravity waves are fast and inertial waves are slow. 
On the other hand, for scales satisfying I773I/I77/I << 1 gravity waves are slow and inertial waves are fast. 
Control of resonances coupling fast waves on small scales and fast waves on large scales can only be achieved 
through a careful analysis of small divisors in resonances (BMN, 1996 a, b, d and Section 6). 

The inertio-gravity wave propagator is the operator solution E(/) (E(0) = Id is the identity) to the linear 
problem (2.3): 


<F(f) = E(t)<F 0 ; $(0) = $ c 


(2.5) 
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where 4>(t) = (U(£), p(t)). The operator E(£) describes propagation of inertio-gravity waves. In the Craya- 
Herring cyclic basis (Riley et al. (1981), Metais and Herring (1989), Godeferd and Cambon (1994) present 
developments in this basis for stratified flows without rotation) linear problem (2.3) restricted on the subspace 
of divergence free vector fields reduces to the following 3x3 matrix for the n-th Fourier component: 

0 0 0 

0 0-1 

0 1 0 

where zero eigenvalue corresponds to the quasi-geost.rophic mode. We use the extended notation 

m ~ [mi, m 2 /a 2 , m 3 /a 3 , 0], 0 < a 2 < 1, 




and similarly for n,k. We introduce the orthonormal basis of the divergence- free subspace for n-th Fourier 


mode: 


P«n = hnTT' F7T.0. 0 ; Pm = I pTFTT. tFFF’ ‘TtFTT' °l i Pm = e 4 = { 0, 0, 0, 1], (2.7) 

L \ n \ \ n I J Ll n l \ n \ N \ n \ mil” I J 

The vectors Pok,POm,Pik etc. are defined similarly. The vectors pok,Pik,p2k are orthonormal cyclic 
vectors for the matrix S restricted on the divergence free Fourier subspace; let P</ be the projection on 
divergence free vectors in the Helmholtz decomposition (for the velocity component): 

|«/| 

— 0) P d^nPln — 4 > nP2ni ^ > d^nP2n — fynPln'i 4*11 — TTT • (2-8) 

|n| 

The pj n are the Craya-Herring basis for the purely stratified problem, already used by Riley et al. (1981). 
In the case / / 0 we use the following orthonormal basis 

QO n — — 1 {N (frnPOn T f£nP2n), Qln ~ Pirn ?2 n — {N 4 ) nP2n f4nP0n) (2-9) 

W(! LO n 

where 

tn = ^,*i = N m+fe n . (2.io) 

The vectors qo n , < / ln and q 2n are orthonormal and form a basis in the space of divergence free vector fields. 
In the case ni — n 2 — 0 (this case corresponds to taking horizontal averages) we have from (2.9) (we have 
used the fact that pi n |n 1 = n 2 =o = (^, 0, 0)) 

Qon = (0, 0, 0, 1); q m = (J=, 0, 0); q 2n = 0,0) (2.11) 

where n = (0, 0, n 3 ) denotes wavenumbers for which ri\ — n 2 = 0. In this paper we use the overbar notation 
for the operation of horizontal averaging 


ni n 3 h 2 h 3 -n{ - h 2 


n 77 / n 


27r 1-21T / cir 


F ^' X ^ 2tt2t r/a 2 J 0 J 0 


F(t, X\, x 2 , x 3 )dx\dx 2 . 
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For the matrix of linear problem 


M„ = iV r S n +/R n (2.13) 

the vectors (2.9) form a cyclic basis since, after projection on divergence-free vector fields via Helmholtz 
decomposition: 


P d Mq 0 n = 0, PdMtfin = -U> n q 2n , PdMq 2n — U> n Qln, (2.14) 

and, as u; 2 = / 2 for n\ = n 2 = 0 

PdMtfon = 0, P d Mq ln = -fq 2 n, PdMq 2fl = fq ln . (2.15) 

Any arbitrary divergence-free vector field can be written as 

Un — Vn ( J0n + V^qin + V~q 2n . , (2.16) 

We shall use the variables V to denote vector of coefficients corresponding to Ujy V„ = [H n ° , Vj , H n 2 ] = 
[Vjj 3 , V(J, = [l 7 , 1 , Vjj 2 ]. We denote by the projection onto q 0n (quasi-geostrophic mode): 

n^ G vHx) = ^i;W n,, ) n 7 ? G ui = 14V.. 

n 

Similarly, we define the projection H AG onto ageostrophic component: 

n J,G U t (ar) = J2(Vn<hn + V^ 2 n)e in ”, n^'uj, = Vjq ln + V^n- 

n 

The action of the linear propagator on the Fourier components E(f) can be written in l 7 - variables 

E(i)[V°, V]„ = exp(« n <M;)[V°, V']„ = [V°, exp(u„iJ)V']. (2.17) 


Here J, Mj, are given by 


We have 


M' = 


0 0 
0 0 
0 1 


0 ^ 

-1 , 

0 J 


J = 



(2.18) 


exp(w n Jf) — cos(w n f)I + sin(cj n /)J. (2.19) 

Obviously, E(/) represents vector rotation in V 1 , F 2 -plane; V° component (called 3DQG) is not affected. 
Note that the relation between U and V variables is given by 

Vn = U> ■ ,o„. Vi = U* • «l„ , Vi = V* ■ « 2 „ ■ (2.20) 
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Clearly, T n 0 * = —V® n for real U(x). 

The Boussinesq equations (2.1) in Fourier representation in V variables can be written in the cyclic 
basis (2.9) as 


dtVJ 3 = -i 


E 


~MM' n v n r 


k-\-m—n,i i ,22 


(2.21) 


where i\,i 2 , i 3 — 0,1,2, M' is the matrix M in C-variables given by (2.18). The coefficients Q^n are 
determined from the equations 


a *1*2*3 

kmn 


{9iik - JW )(?* 2 m ' ?*3n)- 


( 2 . 22 ) 


The coefficients , k + m. = n, are explicitly given by a straightforward computation using (2.9), (2.22). 

For the coefficient describing interaction of 0-modes (quasi-geostrophic) we obtain ( n ' A m' = rqmo — nomi). 

N h' A rhf 

(2.23) 


n ooo _ 

hf kmn 


(■ N 2 ti-m' + f h 3 rh 3 ). 


u k u> m u) n \k\\m\\n\ 

Since skew-symmetric in k,m component of Q° k ^ n makes no contribution to Eq. (2.21), we can use the 
following Q° k ° n ° in in Eq. (2.21): 

N k' Am 7 


^000 
'‘v kmn 


(AC | fh'^ + f 


Noj m \rh\ k' A 


(2.24) 


uj k u> m u) n \k\\fh\\n\ ' ' 07 u) k u n \k\\h\ 

where we have used \m\ 2 u 2 n = A r2 |m'| 2 + / 2 m§. The coefficients Q k °^ n given by (2.24) is the familiar 
coefficient in 3DQG equations written in cyclic basis. Other coefficients in (2.21) have the form 


O 012 - 

V kmn 


fNtj m jn(k' Am 

UJ k U> n \k\\77 


2 


0 Qmkn = 


fN(m 3 k' • m’ - k 3 \m'\ 2 ){h 3 n' • k' - k 3 \hf) 


u k u n \k\\m\\n\\m!\\ri 


(2.25) 


O 021 - 

V kmn 


fNUtnjk' Am 7 

u k u m \k\\m'\\h' 


n 201 

'V mkn 


fNZmZnjk' A m') 2 


(2.26) 




i __ A/ r (C A m' )(h 3 m 3 h l • m 1 + j7? / | 2 |ra / | i ) 


Wjc m’ n' fc m n 


O 101 

i ^vrnkn 


N£ n (k' A m')(rh 3 k' ■ mf — k 3 \m'\ 2 ) 


u) k \k\ m \\m' m/ 


(2.27) 


Q 022 _ N(k' A 7h / )(7V 2 jmf jnf + f 2 m 3 h 3 m' • 7 V 
io k Lo m u> n I A 1 1 m 1 1 n 1 1 m' 1 1 h' I 


0 2°2 pNm 3 (k' A - h 3 h' • k') 

w k w m uin |& [ |m| | n 1 1 m' | |?V | 

In addition, we have coefficients associated with the operation of horizontal 


averaging, where k = (0, 0, k 3 ) : 
n 000 - n?°_° - nooo _ n . n oi 2 _ n 02 i _ n qn __ n 022 _ n 0(U 

^ kmn 'vkm.n M i kmn ~ u > ^ kmn ^ kmn ~ ^ kmn ~ b ? kmn ~ 
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Abo, Q^ 0 = o(WO,i 2 ^o) a nd 


r ) l 02 ,-,201 _ „101 _ „202 

* m. k. n. 


(jOr 


~ t *3 > Qmkn ~ Qmkn ~ Qmkn ~ O ' 


(2.30) 


We note that = V^V^ 1 and then collecting terms in Eqs. (2.21) we define the symmetrized 

/son _ on i nioi /soi 2 _ rjoi2 , n 102 n 021 — n 021 xn 201 n 022 — n 022 n 202 (o q-n 

'vkmn y °6kmn ' Vmfcn> Vfcmn Vfcmn ' q? mkni V kmn — Ka ikmn ' 'vmkni ^ •Ckmn — 'V kmn ' ' a imkn- 

It implies that 


Q o2i D- Ll — D- zz n D- 12 — 

kmn kmn ^ kmn T h/ kmn 


iOII 


*022 

'kmn ^ kmn ^ kmn 


012 

kmn 


Nfim&n r 


*3. 


LO* 


(2.32) 


3 The limit equations describing reduced dynamics 

We write classical rotating Boussinesq equations using the variables V in the basis (2.9); these equations 
have the form (see Eqs. (2.21)) 

d t V n +w„M> b = (B(V, V))„. (3.1) 

Here M 7 is the linear propagator operator corresponding to inertio-gravity waves given in cyclic ^-variables 
by Eqs. (2.17)-(2.18). We introduce this linear propagator directly into nonlinearity in (3.1) using the 
change of variables V = E(-t)v, V n = exp(-o; n M( l t)v n where v = [u°, v 1 ,!) 2 } and M' is defined by (2.18). 
Boussinesq equations (3.1) written in v variables have the form 

d t v — B(t, v, v), B(£, v, v) = E(2)B(E(-t)v, E(-t)v). (3.2) 

Equations (3.2) are explicitly time-dependent with rapidly varying coefficients. The following equations 
describing reduced dynamics are associated with Eqs. (3.2) (BMN, 1996b): 

<9 f w = B(w, w), B{w(t),w(t)) = lim T _> +00 ^~ j B(s,w(t),w{t))ds. (3.3) 

1 Jo 

Clearly, when represented in Fourier modes, the right-hand-side of (3.3) will be determined by resonances 
±w k ±u} m ±tu n = 0 within terms of the type exp(i(±ui k ±u) m ±u; n )s), see (2.17), (2.19). Projecting Eqs. (3.3) 
on the 3DQG mode qo n corresponding to the zero eigenvalue of the linear problem, we find that the limit 
equations (3.3) describing reduced dynamics contain classical 3DQG equations as a completely decoupled 
subsystem, a result already published in BMN (1996b). This confirms the QG decoupling of the reduced 
equations after their projection on the mode corresponding to the zero eigenvalue which was also proven for 
the rotating shallow-water equations and 3D rotating Euler equations (Mahalov and Marcus, 1995; BMN, 
1995; BMN, 1996a; BMN, 1996b). In the case of 3D rotating Euler the reduced equations were projected 
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on the subspace of 2D-3C fields by means of the operation of vertical averaging. Two-dimensional, three 
components refer to fields that have three components and depend on two variables. These are MFI (material 
frame indifferent fields) uneffected by rotation (Ristorcelli et al. , 1995; Speziale, 1989). In the case of rotating 
shallow- water equations reduced equations were projected on the quasi-geostrophic mode and decoupling of 
quasi-geostrophic equations was obtained and published in Mahalov and Marcus (1995) before Embid and 
Majda. (1996). In the important paper by Warn (1986) the same result for shallow-water equations had 
been obtained based on formal two time scales expansions and Ertel’s theorem. To prove the decoupling for 
Boussinesq equations, Embid and Majda (1996) also used ErteFs theorem, which does not give insight into 
the structure of the ageostrophic dynamics in the context of Boussinesq equations. 

Now we ouline a new proof of decoupling of 3DQG equations in (3.3) We refer to BMN ( 1996a, b,c) and 
Mahalov and Marcus (1995) for detailed proofs of splitting in Boussinesq, 3D Euler and rotating shallow- 
water equations. By (3.2) and (3.3) the only active triadic interactions in the reduced equations (3.3) are 
those satisfying the relation ±uik ± ui m ± u n = 0. We note that projection of Eqs. (3.3) on QG mode 
(which corresponds to zero eigenvalue of the linear problem) leads to the additional constraint u> n — 0. Then 
the conditions ±u d= to m ± u) n — 0 and ui n = 0 reduce to 2-wave interactions u>k — u? m . For N 2 ^ / 2 the 
condition cok — u> m is equivalent to the condition |fc / |/|X*| — (m'l/lml (equivalently, <pk — <f>m\ see (2.8)). 
Clearly, reduced equations (3.3) projected on QG mode involve only the coefficients with i 3 = 0 

( n ~ Ar + m). One trivial solution of = u> m is cok = ’m = 0 which corresponds to the QG coefficient 
An important observation is that other terms involving the coefficients Q^ii ( ? 'i / 0 or ? 2 7^ 0) in Eqs. (3.3) 
are annihilated when the resonance condition <j> k = (j) m is used. For the Boussinesq case, the detailed proof 
can be found in BMN (1996b). 

Thus, in the reduced equations (3.3) the total velocity splits into the limit quasi-geostrophic field 
U^(i, x\, x< 2 , x 3 ) satisfying (3.12) and ageostrophic component satisfying equations of the type: 

dtWAG = B2(U*(J), w^g) + B 3 (w^g, w^g)- (3.4) 

In this paper we are interested in the situation when the quasi-geostrophic (balanced) field is present and 
it is dominant. We explore the implications of the existence of the vortical mode (at mesoscales the vortical 
mode represents quasi-geostrophic flows) on the dynamics. The major shortcoming of many existing theories 
is that they neglect the interaction with the vortical (=potential vorticity carrying) mode of motion (Muller 
et al., 1986). The results of the study by Lelong and Riley (1991) provide further evidence that the role 
of the vortical mode in influencing the evolution of strongly stratified flows may be significant and should 
not be neglected. According to Muller et al. (1986), interactions among internal waves and vortical modes 
loom as one of the important questions to be addressed by a strong interaction theory. In BMN (1996 a-e) 
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and here we present such strong interaction theory based on rigorous KAM type mathematical analysis and 
small divisor estimates. 

In BMN (1996b) and in Section 6 it is shown that the main contribution to resonances in the reduced 
equations (3.4) is given by the term B 2 (U^ (t) , w^g)- We arrive at the same conclusion empirically comparing 
orders of magnitude of the term B 2 (U^(Q, w^g) and Bs(w j 4 g,w^g): 

JB, _ (wx £ ,w A q) | w M wei (3 . 5) 

iB 2 (U T (t),w^G)| |U T | 

for synoptic scales at midlatitudes (e.g. Holton, 1992). Then equations become linear in w ag with coefficients 
depending on arbitrary quasi-geostrophic field U \t): 

dt^AG = B 2 (U t (i), w^g)- (3-6) 


The numerical work of Bartello (1996) also shows that resonant 3-wave interactions are of secondary impor- 
tance in the overall picture of interactions when both rotation and stratification are present and Bu — 0(1). 
They do not lead to slow-fast energy exchange and are difficult to resonate. Rigorous proof of this result was 
given in Babin, Mahalov and Nicolaenko (1996b) using small divisor estimates and is discussed in Section 6. 
Geostrophic adjustment takes place via catalytic interaction involving a “rotational” (“vortical”) mode and 
two inertial-gravity modes. As in Bartello (1995), an interaction is “catalytic” in that it does not influence 
slow modes, but serves to transfer fast AG energy downscale. 

The bilinear form B(w,w) in Eqs. (3.3) can be conveniently represented in the cyclic basis (2.9). Now 
we write the resonant (limit) equation for w = w°qo + w 1 q\ -(- w 2 q 2 where w° is simply projection of w 
on QG mode. Then w° can be defined from the quasi-geostrophic equation; w' = (w 1 ,w 2 ) is found from 
inertio-gravity wave limit equation. The quasi-geostrophic equation is given by 


d t w 0 n = Bi(u>°, w°) n , Bi(w°, w°) n 


->■ Qkmn W k™l, 

k+m=n 


G)000 _ 

Vicmn 


Nui m \m\ k' A m' 
u k u) n \k\\h\ 


(3.7) 


and in the case n\ = n 2 = 0 


Bi 


E « 

k+m=7i 


~Wi„w: 


kmn k 


(3.8) 


since Q k ^n — 0 by (2.29). From (2.29), we also prove that p and the horizontal averages of velocities do not 
enter the QG equation (3.7); that there are no or w|f . This is not true for the operator splitting in the 
limit / = 0, Fh « 1, Bu » 1. 

We introduce variables q, U (quasi-geostrophic potential and velocity, not to be confused with the cyclic 
basis vectors q o n , q\ n and </ 2n ): 


q m = w m |m|w^, U ft = N[-k 2 /a 2 , An, 0, 0]q k /uil\k\ 2 . 


(3.9) 
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The quasi-geostrophic stream function is defined by = qkN /<jj1\k\ 2 . Recalling that cj|| k \ 2 /N = \k'\ 2 +/j. 2 kl, 
p = f /N we find a familiar formula which relates i/> and q in physical space 

“( V A + / * 2 ^|)^ = ?• (3.10) 

Using (3.9), Eq. (3.7) is written in the form: 


dj i in — i 'y ^ ( U/j • m) q m . 

k-\-m—n 


(3.11) 


q(t,x) obeys in physical space the 3D quasi-geostrophic equations (Bourgeois and Beale, 1994) 

d 


Dt 


q = o 


(3.12) 


where is the advective derivative, based on the quasi-geostrophic velocity and q is the quasi-geostrophic 
potential. 

Now we describe adiabatic invariants associated with horizontal averaging. Let us denote n as (0,0, 77.3). 
We note that in the case n 1 = n 2 — 0 we have q 0fi = [0, 0, 0, 1] and, therefore, = p n . Now the fact that 
Qkmn — 0 implies -§jPn — 0. Then horizontal averages of p are conserved by the reduced equations. We 
outline the proof of this important new result, that is jj-pn = 0. First notice that Q^mn = 0 from (2.23). 
Also, Qllnn = Qfcmn = Qkmn ~ 0- Then p n , the n-th Fourier coefficient of p{t,x 3) satisfies: 


d_ 

dt Pn ~ 


E 


Q 


120 _ 

kmn 


{w\w 


2 

m 


k-\-m=n ,n ‘ = 0 


U>Wm)\ 


(3.13) 


the right hand side is null on 2- wave interactions <f> /. = <j> m with n' = 0. Clearly, in (3.13) no 3 wave 
resonances are omitted. This result holds for all resonances. The proof is very closely related to the one for 
exact QG splitting (see discussion before Eq. (3.4)) but does not follow immediately from Ertel’s theorem. 
This exact conservation law of the reduced equations correspond to adiabatic invariant of the full Boussinesq 
equations. In particular, it shows that in the asymptotic state (after several periods of oscillations associated 
with wave motions) horizontally averaged buoyancy variable will reach a constant in time value. It confirms 
that horizontal spatial averages of buoyancy variable />(* 3 ) are time independent in the asymptotic state; 
this is frequently done in many investigations addressing the impact of vertical variability of buoyancy on the 
dynamics (e.g. Howard, 1972; McWilliams and Weiss, 1996; Doering and Constantin, 1996). Below we find 
the impact of buoyancy phases associated with vertical shearing of this invariant on ageostrophic dynamics 
in non-hydrostatic situation. Similarly, we prove that U 1 , U 2 undergo rigid D-rotation in the Bu = 0(1) 
limit (this result is true for all resonances including 3 wave resonances) and they are true adiabatic invariants 
only if Q — 0. 
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The limit equations for = [w^w^] include w j n , and already found w®.. Therefore, taking into 

account k , m -symmetry, we always label w° variable by k; therefore to write the resonant equations it suffices 
to consider , Qfcmn > (recall that all 1-wave and 3-wave interactions are non-resonant, cf. 

Section 6). The resonant relations reduce to u> n = uj m ,n = k + m for the “symmetrized” Q % kmn defined 
in (2.31). The original Boussinesq equations for v' = fl^^v can be written in the form 


d,v' n = -iexp(ui„iJ) V v° k Q kmn exp(-u> m tJ)v' m +Il£ G Q + (Nt,v) n . 


(3.14) 


Here the matrix Q kmn is defined by 




Q on 


Q 


kmn 

012 

kmn 


Q021 


Q 


kmn 

022 

kmn 


(3.15) 


Qfcmn e ^ c - are gi ven above in (2.31)-(2.32), J is given by (2.18); Q (Nt,v) n are non-resonant terms. Taking 
into account (2.19), we see that resonant interactions are generated by squares of sines and cosines and are 
given by the matrix Q = (Q — JQJ)/2, namely: 


Q kr 


cion , <9022 


_n012 , Q021 

kmn ' '‘v kmn 


(3.16) 


\ _/O UJ1 n 012 nun i /qi 
\ '“v kmn ' kmn ’ ^ kmn ' *‘v kmn ) 

Computing explicitly resonant terms in (3.14) we obtain reduced resonant equations for the limit w, with 


w„ = iw„ , w 


dt.w. 


E 




k 3 pk Jw m - i 


E 


'2(0 Qt 


w 

mn y m 


(3.17) 


k -f- m ~ n , k* = 0 


A- + m = n , fc ' 0 


|n| 


and 


We recall that the resonance condition co m = u> n is equivalent to the condition 

= ( or l£m| = |^n|, see (2.10)). In Eqs. (3.17) w°(t) is an arbitrary solution of quasi-geostrophic 
equation with initial data projected on QG fields and Q kmn are geometric coefficients given above. 


In Eqs. (3.17) 


N4>m <t>r 


k 3 p k represents phase correction impact of vertical shearing of p k on ageostrophic 


dynamics. Clearly, = ik 3 p k . This phase correction implicit in (3.17) is equivalent to the spectral 

Doppler phase shift. It can be calculated explicitly as a solution of the linear problem: 

N (f) m (f) n 


d t w r 


E 


-k 3 pk Jw r 


(3.18) 


In particular, in (3.18) for all a 3 the condition <f) m — <f> n is equivalent to m 3 = —713, m' = n' (for m = n, 
Pk = 0). It implies that k 3 — 2n 3 , <p m .<j) n — <j>^ in (3.18). With the above remark on the convolution in 
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(3.19) 



Here 2n 3 p (o,o, 2 n 3 ) is vertical shearing of horizontal spatial average of p which is an adiabatic invariant (time 
independent) in the asymptotic limit. 

The term containing vertical shearing of p k in (3.17) is associated with the coefficient given 

by (2.32). It involves the vertical shearing operator d/dx 3 (multiplication by* k 3 in Fourier space). In 
this respect Eqs. (3.17) should be useful in detailed investigations of vertical mixing in times of high shear. 
Doppler effects such as shifting the frequency of an internal wave can be measured experimentally through 
Eqs. (3.21). We have shown the existence of a. statistical dephasing effect, induced by turbulent processes on 
inertial-gravity waves. Turbulence acts to renormalize both frequency and viscosity of waves (Legras, 1980; 
Carnevale and Martin, 1982). Legras (1980) has investigated turbulent phase shifts of Rossby waves coupled 
with QG turbulence. Using EDQNM theory and numerical simulations he showed that in the turbulent 
domain turbulent frequency shifts can be as large as 30 % of the wave frequency given by linear theory. We 
effectively generalize this work to the most general interactions of inertio-gravity waves with 3DQG turbu- 
lence. In BMN (1996a) we presented formulas for frequency shifts induced by turbulence on inertial waves in 
the context of 3D rotating Euler equations. The turbulent shifts for rotating shallow-water equations with 
the (3 term were analyzed in BMN (1996 c, e). 

The equations (3.17) for w^g are also invariant under vertical averaging, k 3 = rn 3 = n 3 = 0; this follows 
from the equivalent condition for two wave resonances, L 0 m = ui n implies |m 3 |/|m| = |n 3 |/|n|. Denote by 
w AG the vertical average of w AG . Careful inspection of the coefficients in (3.16) shows that they all 

reduce to zero, except for: 


Of 1 - 1 - = Q W- 

kmn kmn 


k' A m! 

1*1 


(3.22) 


Moreover the Craya-Herring basis reduces to q Q = = p Q =, q^= = (0,0, -1,0), q 2 ~ = (0,0,0, 1). Hence W AG 
reduces to the vertical averages of vertical velocity and buoyancy variables. Vertical averaging is usually 
denoted as the barotropic component (Bartello, 1995). The above establishes that these quantities are 
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purely passively advected by the 2D-2C vertically averaged quasigeostrophic velocities. The linkage with the 
dynamical Taylor-Proudman theorem at Bu — 0 is remarkable, as U3 and ~p are indeed passively advected 
by the 2D-2C vertically averaged Euler fields in that context too. It confirms that barotropic components 
ii 3 and p of the field take no part in the geostrophic adjustment process. 

Further, for almost all Bu and domain aspect ratios a 2 , B 2 (wqg', w^g) in (3.17) splits in Fourier space 
into uncoupled, restricted interaction operators on 4 rays families 


^ mi ^ 


H- 

0 

0 


n, \ 

m 2 

= A 

0 ±1 0 


«2 

V ms ) 


v 0 0 ±1 j 


( " 3 ) 


(3.23) 


for n , 3 / 0. This is obtained by further reducing resonances ui m — u n . In Eqs. (3.17) together with 
the condition (3.23) direct cascades of energy are allowed for w ag through B2(wqg,w^g). Wave energy 
cascades toward smaller scales and is subject to strong dissipation in AG mode. Vertical shearing operations 
in B 2 (wqg ]W^c) are conveniently expressed in terms of divergent velocity potential X . We recall that 
X and -0 are related to U h by the formula U h = e 3 x V h^P + V hX- The form B 2 (wqg, w^g) takes an 
especially simple form in the Craya-Herring cyclic basis. In this basis the first component of wag is simply 
the divergent velocity potential X — ( — A)^ 1 d Q~f - or in Fourier space 


Xn — 



;WZ 


(3.24) 


where h' — («i,n 2 /a 2 ), and the components of the full field in the Craya-Herring basis are 

Here w® is the component of the quasi-geostrophic mode and tn* , w~ are the ageostrophic components. 

The coefficients <5j.“ n = Q^ n , Qf„\ n in (3.15)-(3.16) explicitly yield terms Xm = “* 

Xn — which reflect vertical shearing dynamics, and contribute to imbalance dynamics. These 

terms come from Q ( ^ n , Qf^ n , QmL S iven b y Ec l s - (2.25)-(2.27) and are carried over to the 

Qkmn > Q^mn hi Ec l s - (2-31) and in Eqs. ( 3 . 1 5) - ( 3 . 1 6 ) . The important role of the divergent velocity 

potential \ in the AG reduced Eqs. (3.17) can be compared with that in the classical balanced models (Gent 
and McWilliams, 1983a, b). There a classical expansion in small parameters is carried. In some sense, our 
QG-f-AG Eqs. (3.17) describe near balanced and some unbalanced regimes. 

Wave energy cascades towards smaller scales along rays in Fourier space and is subject to strong dissipa- 
tion in the ageostrophic mode. We have an upscale energy cascade via (0, 0, 0) triads corresponding to 3DQG, 
downscale cascade of energy via (0, i 2 , is), i 2 7^ 0, i 3 / 0 triads corresponding to catalytic interactions. Tri- 
ads (0,7.0,23), ?‘o 7b 0, i 3 7^ 0 flush the gravity wave energy downscale in a nonlinear geostrophic adjustment 
process. When simulations are initially balanced, gravity modes act to damp large scale rotational modes 
via transfer into intermediate scale gravity modes, and via a subsequent downscale wave cascade involving 
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the catalytic interaction. Effective eddy viscosity on fast modes is larger than that on slow modes confirming 
Bartello et al. (1996). Spectral eddy viscosities for QG and AG fields with explicit dependence on the rota- 
tion/stratification parameters are calculated. They correspond to 2- wave resonances on the linear diffusion 
operator. Let and be the kinematic viscosity and the heat conductivity, respectively; the ratio vxfv 2 
is known as the Prandtl number. Then effective spectral eddy viscosities uqg and vag are given by 


VQo{n) = " 2 + ( "‘ ~ "Vp+Ur VAaln) = 1/1 + ( " 2 " ^ G'p+Ug 


(3.25) 


where p — f/N. Let’s consider the case of a large Prandtl number vi/v 2 >> 1- Then Eqs. (3.25) can be 
approximated by the following expressions 


"« o(B) = VAa(n) = ^ FF+kr 


(3.26) 


Eqs. (3.26) clearly show that QG and AG modes are effected by viscosity differently. For example, vag{k) 
increases if vertical shearing increases (large ns) while vqg(ti) decreases. The dependence of effective eddy 
viscosities on p = f/N is also of interest. In the stratification dominated regime (p — » 0) we have UQa{n) -> 
J'l, ^ag(w) — > 0. The fact that oag{^) — > 0 is reflected in lack of control of AG dynamics (vertical shearing) in 
the strongly stratified limit. The situation changes dramatically in the rotation dominated regime p — » + 00 . 
In this situation we have uag{N) — > tq, UQc(n) — > 0. The fact that i>ag{N) approaches a nonzero finite value 
allows us to control AG dynamics in this limit. Remarkably, UQa(n) — >• 0 does not spoil the situation since 
QG dynamics is globally controlled even for zero viscosity (Bourgeois and Beale, 1994). In the physical case 
Pr = vi/v 2 1 and naively adding the usual viscosity to 3DQG equations is incorrect in the geophysical 
limit. 


As noticed by Metais and Herring (1987) and Metais et al. (1996), the presence of stratification yields 
the formation of very strong vertical variability which tends to destroy the vortices vertical coherence. In 
stably stratified, rotating flows, quasi-two-dimensional organized vortices structures are only observed in the 
presence of a well-defined horizontal density (temperature) front leading to baroclinic instability. 

The impact of even small viscosity onto ageostrophic dynamics in (3.17) is strikingly different from that 
on the quasi-geostrophic equations (3.12) confirming the predictions of Metais and Herring (1989). In (3.12) 
modified with vqg dissipative small scales interact nonlinearly with an anisotropic inertial range with inverse 
cascades and intermittencies in the quasi-geostrophic turbulence. Inviscid ageostrophic dynamics are driven 
by the chaotic quasi-geostrophic field, albeit on restricted rays in the Fourier wavenumber space (but no 
restriction to the direct energy cascade). Energy is conserved separately on each family, but not enstrophy. 
Direct cascade of energy is enhanced along each ray family, in contrast to the inverse cascade for the QG 
component spectrum. This is confirmed by numerical observations of Metais et al. (1996). 
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4 Intermediate asymptotic regime of strong stratification and weak 
rotation 


In this section only £ 2/02 will be denoted by &2 and k^/a^ by k 3, |A‘j 2 = kf + k^/ao + k^/a^, \k'\ 2 = k\-\-k 2 / a?>. 
Intermediate asymptotic regime of strong stratification and weak rotation is conveniently characterized in the 
Craya-Herring cyclic basis using expansions in a small parameter p = f/N = Jq/Nq. Equations describing 
balanced dynamics with control of vertical shearing are obtained by expanding Q^mn = ( k ' A mf) 

in Eqs. (2.24) in powers of p. We have 

“*1*1 = jN*\k>\* + Pkl = N^\k’\i + fki = AT|*'|(1 + + 0(P)), (4.1) 


N _ 1 


— (1 

| k'\ { 2|Ar'| 


r /i 2 + 0(p 4 )) 


(4.2) 


with similar expressions for w m |m| and to n \n\. Here 0(p 4 ) denote terms in Taylor series of order p 4 and 
higher. Expansions (4.1) are valid provided that « 1. We note that l v ~ 1 / 1^3 1 and lh = 1/|&'| are 

the vertical and horizontal spectral length scales, respectively. Then the condition jp^p 2 << 1 is equivalent 
to the condition Bu s = N^l 2 / JqI^ » 1 . Here Bu s is the spectral Burger number. 

Substituting expansions (4.1) in the expression for Q^mn we obtain after cancelling N 2 and suppressing 
terms of order p 4 and higher 


QZnU 0 = 


1 * 11 ’ 


Ik 1 A 




,/|2 


1*1 




Clearly, 


(4.3) 


«„ ( °) = ||p| ( *' A "')- (4.4) 

The appearance of the operator kl/\k'\ 2 etc. in Eqs. (4.3) is not surprising since it is related to the inversion 
of (3.10). Indeed, for small p we have 

d 2 d 2 

* = -(Vi + r 1 ? = -V,- 2 (l - / /2 V/T 2 ^|) + 0(p 4 ) (4.5) 

where the operator — is —k^/\k'\ 2 in Fourier space. 

Eqs. (4.3) show that the regularizing effect of rotation involves the inverses of horizontal Laplacian. 
The effect of these operators on the dynamics is to decrease the horizontal scale. This is in fact what is 
observed in numerical simulations investigating the effects of rotation on turbulence (McWilliams, private 
communication) . 
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Eqs. (3.7) now take the form 


dtW° n = -i J2 Q°kmn{M)w 0 k w° m '(4.6) 

k+m=n 

where Q k ^ n (fi) is given by Eqs. (4.3). For every /i ^ 0 the coefficients Q^mnif 1 ) h ave explicit dependence 
on vertical variability parameters {k^, 1113,113 in Fourier space). For g = 0 dependence on vertical variable 
disappears and we obtain complete decoupling of dynamics at every vertical layer with the vertical coordinate 
X 3 appearing only as a parameter. At any given vertical level the dynamics is described by 2D Euler 
equations with different pressures at every vertical level. There is no control of vertical shearing which can 
grow exponentially. We obtain quasi-2D unbalanced dynamics of vertically stacked pancakes which was first 
described by Lilly (1983). For // f 0, as shown below, Eqs. (4.3), (4.6) have a conservation law which allows 
to control vertical shearing for all times. Our results reveal the regularizing effect of weak rotation which 
introduces vertical coherence; the mechanism of the regularization is through coupling of dynamics in the 
vertical direction. This regularization differs from regularizations based on vertical viscosity. 

We note that the coefficients Q k ( ^ n (g) have the following important skew-symmetric property 

QlZ(v) = -QnZ(h). (4.7) 

Eqs. (4.6) have a. conservation law which simply expresses the fact that potential vorticity is conserved if 
written in q variables. Multiplying (4.6) by w Q _ n and summing over n we obtain 

ftE>“WI 2 = 0 (4.8) 

n 

since 


-* J2 Qkmn{fi)u>kWmW-n = J2 QTmn{l*) w k w m. w n = 0 - ( 4 - 9 ) 

k^m—n k-\-m~\-n = 0 

In (4.9) we used the fact that Q k ^ n (ft) is anti-symmetric (Eqs. (4.7)) and is symmetric under 

interchange of k and n. The conservation law for w® (t.) implies conservation of energy for the projection of 
the field on 

2 

9on{n) = (1 - ^p/l 2 )Pon + /“j”7|P2n + 0{fl 3 ). (4.10) 

Then using the conservation law (4.8) and following Bourgeois and Beale (1994) one can prove global existence 
of solutions of Eqs. (4.3), (4.6) with full control of vertical shearing for all times. 

The coefficients (p) given by Eqs. (2.25)-(2.28) can be expanded in powers of fi (see Appendix). 

Then we obtain 


QllKh) = Qt 


000 

rnn 


( 0 ) + 


\m 


w\ 


(k' An') (- 


H 


m' 


l*'l 


) + 0(g / 


(4.11) 
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A Oil / x _ Aon /M 1 2 k li k ' A rn')(n 3 m 3 n' • m! + |m'! 2 K| 2 + n 3 m. 3 k' ■ m! - n 3 k 3 \m '\ 2 ) 

Gkmn ' c efcmnV'-V r> V 1 




, (4.12) 


A022 / x _ A022 /n\ _ 1 2 k ' A m ' ( k l , m 3 , n l x , 

yimnl J 2^ |*'| V'! 2 + K| 2 + |«'| 2 + 

2 (&' A m!){m 3 n 3 m. 1 ■ n' + m 3 k 3 \n'\ 2 — m 3 n 3 n ’ ■ k') 




/,l m'||n / ||m / ||n / | 


+ 0(n ), 


(4.13) 


A 012 / x _ o rn 3 n 3 (k' A m') 2 + (m 3 k' ■ m! - k 3 \m'\ 2 ){n 3 n' ■ k' - k 3 \n'\ 2 ) t ^ 

Vfemn KN ~ ~*L l 1 7 /ij /II fl /|| /| r U(fJ j, 

| A?' | j r?' 1 1 y?2 1 1 m' 1 1 « | 


(4.14) 


A021 - 1„ m 3”3( y Am ') 2 , nln 3s 

Q km M-4 ||fc , +0(n ) 


where 


(4.15) 


QZ°n( 0) = |4jj4|(*' An'), QfA(O) = <S(0) = QK„(0) = 0, 


(4.16) 


Q*m»(0) = 


(A;' A m')(n 3 m 3 n' ■ in' 4- Im'pln'l" + n 3 m 3 k' • m! — n 3 k 3 \m 
/ 1 1 / 1 1 / 1 


/| 2 ) 


(4.17) 


Our regularized system for w n = (w^w' n ) t w' n = (w„ , w 2 ) follows from Eqs. (3.7), (3.16)-(3.17). It has the 


form 


d t u, n = “* Y1 Qk°mn(^) w k w m^O{f.l 4 ), 


k+m=n 


(4.18) 


<9*w, 


* hp k Jw m - i ]T Wfc(<) Q*mn(^)Wm +0(jU 4 ). (4.19) 


4>m = 
fc -f m : 


*£m — ■ 4>n 
k TYl — 11 


where 


~ / ) _ ( Dkmn(v), Gkmn(p) |_ n / \ [ ^ M i e / x [ ® ^ 

Q kmnif^) — I I — D kmn\[ ^) I I + O-mn (/^) I 

\ -Gkmn(v), D kmn (n ) / \ 0 1 / \ -1 0 


(4.20) 


Dkmnin) — Qkinn (l J ) + Qkmnil 1 ) ■ G kmn (f() — Qkmri (/0 + Qkrnnik) 


012 


(4.21) 


with components of the matrix Q fcmn (a 4 ) given by formulas above. For every fixed value of the parameter p 
vertical shearing is controlled in the reduced equations (4.18)-(4.20) for all times. The system of equations 
(4. 18)- (4.20) will be called the regularized Boussinesq system in the regime of strong stratification and weak 
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rotation. It allows us to study the regularizing effect of weak rotation on quasi-2D dynamics of decoupled 
pancakes and its impact on AG dynamics. For /r = 0 we obtain 

G*m,(0) = Ofcmn(O) - <&(«) = 0, D kmn ( 0 ) = £?“”„( 0 ) + <?““„(()) (4.22) 

where and Qkmn (0) are given in (4. 16)- (4. 17) . We recall that nd and re, 2 , are projections of the 

total field on the divergent velocity potential (x) and the geostrophic departure mode, respectively. The 
matrix Q fcwn (0) is diagonal and equations (4.19) for w * and w 2 are coupled for p = 0 only through the 
phase term /? 3 pk Jw m which is associated with vertical shearing of horizontally averaged buoyancy 

variable (adiabatic invariant). For fi ^ 0 the modes w * and w 2 are coupled as can be seen from Eqs. (4.20). 
Non-diagonal terms are proportional to Gkmn{p)\ Taylor series expansion of Gkmn{p ) starts with terms of 
order //. 

In order to obtain a simple reduced system suitable for practical numerical implementation, terms of 
order //. in Eqs. (4.19) are neglected. Then two scalar equations for w 1 and w 2 - in (4.19) are coupled only 
through the phase term and can be solved in parallel. 

A simplified reduced system consists of three prognostic equations: 

dtWn = -i Qkmn(v)™° k W? n , (4.23) 

k+m=:n 


d t w' n = i 




hpk Jw m - i «$(*) D kmn {0)w m . 


(4.24) 


•Pm — <Pr 

k + m = 


<P ni = <Pn 
k + m = n 


The scalar coefficient D k mn{ 0 ) in (4.24) is given by 

A022 


D kmn ( 0) = Ojj“ (0) + <&( o) = 2 ^ ' >»' + MVP 


<j) n . The first term k in (4.25) represents advection 


In the derivation of (4.25) we used the relation <p t 
in physical space. Now we analyze the role of the second 0-order term 

_ nzm^n' ■ m! + |j 7 i'| 2 |?r / | 2 _ In' A ra '| 2 + (n ■ m)(n' • m!) 


(4.25) 


m) 


\m' 77 .' m 72 


j 777 1 1 77. 1 1 777. 1 1 77 1 

The coefficient. 7 ( 77 , 777 ) can be easily computed for 4 ray families. For example, we obtain from (4.26) 


(4.26) 


7 ( 77 , 777 .) — 1 if mi = A? 7 i, 7 ?7 2 = A 77 . 2 , m 3 = A 77 . 3 ; 


(4.27) 


7(? 


777 2 + 777^ — 777,3 


if mi = At 7 i, 777 o — A77 2 , m3 = — A77.3. 


(4.28) 
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The appearance of a wave type operator — ; in Eqs. (4.24) is due to polarization induced by waves; 
it corresponds to wave dynamics on the reduced family of 4 rays. 

The advantage of prognostic equations (4.23)- (4.24) over the diagnostic vertical velocity (omega) equa- 
tion (1.17) is in improved accuracy. Near balance and some unbalanced regimes can be treated by Eqs. (4.23)- 
(4.24). These regimes cannot be treated by the diagnostic equation (1.17) where filtered initial data is 
required in order for (1.17) to be a reasonable approximation of (1.16). 

For n = 0 Eqs. (4.18) coincide with the familiar quasi-2D Euler system which can be seen by introducing 
variables q and U q 2 d (quasi2D potential and velocity): 

<lm — Uq 2 £>,& = [~ & 2 , & 1 , 0, 0\qk/\k'\ 2 . (4.29) 

In this notation (3.7) is written in the form of 2D Euler equation which depends on X 3 as on a parameter: 

dt(jn. — i ^ ^ (Uq 2£>,A- m) (jm ■ (4.30) 

k-{-m=n 

Then the velocity UQ 2 i>(f , # 1 , * 2 , X 3 ) satisfies quasi-2D Euler system 

dt^Q 2 D + (Uq2D + Uq2d) ' V / Uq 2 J[) = —Vp, V' • Vqod = 0 (4-31) 

which depend on x 3 as a parameter, V' = [<9 i,< 9 2 ]. In Eqs. (4.31) U q 2 d denotes horizontally averaged 
velocity which is an adiabatic invariant of 3D Boussinesq equations in the strongly stratified limit (in the 
absence of rotation). 

5 Phenomenological analysis at asymptotic limit of strong rota- 
tion/stratification in the Burger number of order one situation 

At. asymptotic limit of strong rotation/stratification, the existence of two disparate time scales indicates a 
phenomenological analysis similar to that of rotation (Zhou, 1995; Mahalov and Zhou, 1996) may be appro- 
priate. The aim of this approach is to estimate the averaged effect of rotation/stratification on turbulent 
energy transfer. The introduction of the anisotropic time scale based on the aspect ratio parameter in the en- 
ergy spectrum is an improvement over our previous phenomenological analysis of rotating turbulence. In the 
context of the quasi-geostrophic equations for a Boussinesq fluid in a uniformly rotating and stably stratified 
environment, McWilliams et al. (1994) showed that their solutions exhibit significant anisotropy associated 
with the emergence of many long lived coherent vortices that control the flow evolution. Anisotropy of 
quasi-geostrophic field impacts on the ageostrophic gravity wave field through Eqs. (3.17). 
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Among more mundane immediate consequences, exact operator splitting of the reduced equations (3.3) 
reach the very roots of the mechanisms of wave-vortex interactions. Let us designate by ( ) the ensemble 
averaging for any field TJd, and by = Ut — (U^) the fluctuations. Then the Reynolds stress operator 
becomes: 

(WqG'F, Wqg,f) + <W ag,f,W ag,f) + 2(Wqg iJ p, Wag,f)- (5.1) 

The last tensor we shall designate as the anisotropic phase coherence tensor. It correlates the fluctuations of 
the ageostrophic W^G-held with the quasi-geostrophic WQG'-field. The anisotropic phase coherence tensor 
is the key player in the control of rapid 3-D pressure fluctuations. The dependence of the full Reynolds stress 
tensor on the intrinsic mean vorticity does not vanish in the limit of strongly rotating/stratified turbulence, 
as neither the field G, nor the “anisotropic phase coherence tensor” vanish. Long-lived phase coherence 
is an important part of turbulence (Bartello and Holloway, 1991; Herring and McWilliams, 1985). 

In order to infer the form of the inertial-range spectrum, it is necessary ta estimate the magnitude for 
the triple correlations. In general, 73, the time scale for decay of triple correlations which is responsible for 
inducing turbulent spectral transfer, may depend on any relevant turbulence parameters. Because energy is 
conserved by the nonlinear interaction and a local cascade has been assumed, c is independent of k. Local 
cascade also implies that e is explicitly proportional to r 3 and depends on the wavenumber and on the power 
of the omni-directional energy spectrum. A simple dimensional analysis leads to 

e - A 2 r 3 (k)k 4 E 2 (k) (5.2) 

where A is a constant. We recall that u = yj N{j . As in Section 4, we defined the vertical and 
the horizontal spectral scales as l v — p^, l) x ~ jpy. If 77 = l v /lh — l^l/l&sl is the ratio of these length scales 

then “* = \/ N o^Ti+folAi- 

In the strongly rotating/stratified case when both effects are of the same order (Burger number of order 
one situation), the time scale for r 3 (k), the decorrelation of the triple velocity product, is the controlling 
parameter to influence the energy transfer process. In a regime of high Reynolds numbers and low Rossby 
and Froude numbers, turbulence is characterized by two disparate time-scales: a short anisotropic time scale 
associated with the rotation/stratification frequency t^n = 1 / yj A' ,y + /q (/o = 2O 0 is the Coriolis 
parameter) and a nonlinear time scale. We find that a direct application of r 3 = r^iv resulted the energy 
spectrum for turbulence subject to strong rotation/stratification: 

E(k) = Cnjv(Gi e) l/2 k- 2 . (5.3) 

We now consider the number of the non-dimensional parameters needed. For the turbulence in equilibrium 
the Rossby/Froude numbers are the only relevant parameter controlling the effects of rotation /stratification 


30 



on the flow. However, for the non-equilibrium situation, anew non-dimensional parameter such as t) is 
required. Here the dimensionality is considered by introducing the aspect ratio. Our analysis suggests that 
the energy transfer process in the limit of (r^ t) —> oo and small spectral Rossby/Froude number (strong 
rotation/stratification) is as follows. There is inverse energy transfer by 3D QG component (McWilliams, 
Weiss and Yavneh, 1994). In the meantime, there is also a direct energy cascade governed by the equation 
for the inertio-gravitv wave component. The energy spectrum of the full field Ut is given by E(k) = 

Following the usual assumption of EDQNM, we consider that the lifetime of triple correlations in rotat- 
ing/stratified turbulence might be more accurately treated by taking into account the possibility that these 
correlations decay because of the influence of both wave propagation and nonlinear interactions. The simple 
choice 


r 3 (k) r n i(k ) + mN (k) (5 ' 4) 

satisfies the appropriate limiting cases: 73 ( fc ) — » r n i without rotation/stratification and 73 (Ar) t^n with 
strong rotation/stratification. The introduction of the anisotropic time scale based on the aspect ratio pa- 
rameter in the energy spectrum is an improvement over our previous phenomenological analysis (Zhou, 1995; 
Mahalov and Zhou, 1996) since now the model can distinguish the anisotropic nature of rotating/stratified 
flows. The generalized inertial range energy spectrum is 


E(k) = Z ^-1/3 f 2/3 / --5/3 j 


(5.5) 


where Z is given by 

z =\ U/Y + \l-Y + 2\/y 2 + 4Zo ) , 


(5.6) 


Y = 






(5.7) 


The parameters A = C K 3 ^ 4 , k cin = and Zq = A 2 / 3 r^(efc 2 )~ 1//3 = [ ^p^ ] 2 / 3 . The strong 

rotation/stratification limit then leads to Cu n — 1/vT = 1.22 — 1.87 for the typical range of Kolmogorov 
constant. These equations reduce to the classical Kolmogorov “-5/3” spectrum when Z o — >■ 0 (so that Z — ¥ 1), 
and to our rotation/stratification modified “-2” spectrum when Zq » 1 (so that Y — >• 0, Z — * Z\j A ). For 
intermediate rotation/stratification rates the spectrum varies smoothly between these two limiting forms, 
according to the increase of the controlling parameter Zq with increasing ratio kQ,]^/k. This confirms the 
numerical simulations of rotating stratified turbulence by Metais et ctl. (1996). In their numerical simulations 
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the geostrophic energy spectrum was steeper than k~ 3 at the small scales, whereas the ageostrophic energy 
spectrum was much shallower. As a result, the large scales are dominated by geostrophic energy, while the 
reverse is true of the small scales. There is a spectral gap between QG and AG spectrum with AG spectrum 
being shallower than the typical k~ 3 QG spectrum; it varies smoothly between k~ 2 and k ~ 5 / 3 which is in 
agreement with numerical simulations of Metais et al. (1996), Bartello et al. (1996). 

The equation which determines the spectral eddy viscosity, ^(A 1 ), is given by 

e — i/T(k)S 2 (k) = v T (k)E{k)k 3 . (5-8) 


Therefore, the rotation/stratification dependent eddy viscosity can be estimated as 


v T {k) = u 0 [ 


m,i,2 1 


1 + y ~ 1/2 


(5.9) 


where is a constant and y = k 3 E(k) / N 2 . The inertial range wavenumber k can be related to the turbulent 
kinetic energy K and dissipation rate e 


k = 


E{k) 

R3/2 


(3Cfc/2) 3 / 3 . 


For the inertial range spectrum, Eq. (5.9) can be rewritten in physical space: 

, . , K 2 1 


UT[X 


1 -F 0.36A7(r^]ve) 


(5.10) 


The eddy viscosity above is only appropriate for homogeneous turbulence without mean velocity gradients. 

The spectral time scale is an important measurement. The Kolmogorov hypothesis implies that the 
energy-containing range excitation does not affect energy transfer within the inertial range. Therefore, the 
average rate of energy dissipation is identified with the rate of spectral energy transfer and the rate of energy 
input. The nonlinear (or eddy turnover) time-scale, r nl (k) = [k 3 E (k)]~ 1 / 2 } j s then equivalent to the spectral 
transfer time, r s . Zhou (1995) has shown that 


r s {k) = 


1 I Tm{k)] 2 

A 2 nusr(k) 


(5.11) 


or more generally, 


r s (k) ~ 


A 2 T 3 (k) • 


(5.12) 


Therefore, the time for the spectral energy transfer is increased to a value greater than r nt ; thus the nonlinear 
energy transfer is suppressed by rotation/stratification. 
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6 The devil’s staircase of convergence results 


In this section we give a short presentation of new regularity and strong convergence results for rotating 
Boussinesq equations in the Bu = 0(1) regime; these are substantial improvements over the corresponding 
ones in BMN (1996b) in that much less regularity is required for the initial data, with the Sobolev space 
H 9 being the worst and // 4 being the best case. Here H a designates the space of functions with square 
integrable derivatives up to order a. Any strong convergence result uniform in the parameters Ro a or 
1/N cannot circumvent control of small divisors and sharp estimates of both near 3- wave and near 2- 
wave resonances. The o(l) non-uniform convergence results in H 3 of Embid and Majda (1996) cannot be 
extended to 0(1 /N) uniform error without control of such near-resonances. Here we investigate the density 
and probability of both 3-wave and 2-wave resonances as function of the three geometric parameters: 

Bu = N 2 al/f, 0 2 = l/a 2 2 , 9 3 = 1 fa\. (6.1) 

For rare, non-generic values of these parameters, the limit equations for the ageostrophic field wag can be 
nonlinear: 


dt^AG = B 2 (w QG (t),WAG) + B 3 (w AG , w^g); (6.2) 

the bilinear form B 3 can easily be computed in the Craya-Herring cyclic basis and is not detailed in this 
section for conciseness sake. In B 3 , the domain of summation K* , (k,m) E K* , k + m = n, depends on 
(Bu, 1/a 2 , 1/aff) = (Bu,9 2 ,9 3 )] that is K * = K* (Bu, 9 2 , 9 3 ). For every fixed 9 2 = 1 /a 2 and 9 3 = l/a, 3 , 
the summation set K*(Bu, 9 2 , 9 3 ) is not empty when Bu E @^(^ 2 ,^ 3 ); the set 0^ is very thin, namely 
it is countable. We call it a strict 3-wa.ve resonant set. In Corollaries 6.1 and 6.2 we give an estimate of 
the very small probabilistic measure of near 3-wave resonances. We mostly study the typical case when 
3-wave interactions are absent, that is Bu (f Q*b(9 2 ,9 3 ) and B 3 (wag', wag) is identically zero. When 
Bu E 0]g (9 2 1 9 3 ) is strictly resonant, B 3 is non-zero and strongly depends on Bur, the sets K* (Bu,9 2 ,9 3 ) 
with different Bu do not intersect (highly nontrivial result from our study of the small divisor problem). 
This implies that the operator B 3 depends on resonant Bu dis continuously, every point Bu E 0jg(#2;$3) is 
a point of discontinuity of the operator B 3 . Since B 3 is not zero, solutions of the limit system with general 
initial data discontinuously depend on Bu as well. As solutions of the original rotating Euler-Boussinesq 
equations depend on Bu continuously (on a small time interval [0,Ti]), the convergence to solutions of the 
limit equations cannot be uniform in Bu, a 2 , a 3 . This paradox we call the “Devil’s staircase of convergence 
results” . 

For Bu 0^, Eqs. (6.2) reduces to the linear equations for catalytic interactions (3.16)-(3.17). At the 
same time, the equation strongly depends on 9 2 = 1/a 2 , since the 2-wave resonance condition u> m = u) n 
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is equivalent to jra'I/lm^ = |wVl n 3 l- F° r every irrational 9 3 , Bo splits in Fourier space into uncoupled, 
restricted interaction operators on 4 ray families as in (3.23). For every resonant rational point 62 G @ 2 , B 2 
includes more interactions (much larger, but finite number of Fourier rays), so indeed rational 62 are also 
points of discontinuity for B 2 and further contribute to the Devil’s staircase of convergence. An important 
observation is that % wave resonances are controlled by Q 2 only, not by Bu or a 3 . This follows from the fact 
that u) m = u n implies | m 3 1 / 1 m' \ =■ |« 3 |/|h/|. In contrast, 3-wave resonances are controlled by Bu uniformly 
in $2 and 9 3 \ although strict resonant values Bu = Q* B ( 92 , 9 3 ) depend on $2 and 9 3 , we prove that the 
estimate of the measure of almost resonant Bu does not depend on 9 3 and is uniform in # 3 . The sensitive 
dependence of convergence results on the parameters Bu and 02 , a 3 was missed by Embid and Majda. In 
fact: 

Lemma 6.1: 

Let N 2 — f 2 (that is Bu — a 3 ). Then there exist no 3- wave resonances for 1/9 < a 3 < 9. 

The main convergence result of this section (Theorem 6.1) shows that the convergence is uniform and 
the error is of order (l/p 3 + 1 /^ 3 )/^ when 9 2 6/ Of, Bu (f O 5 3 (^ 2 ^ 3 ), with the Lebesgue measure 
meas(Qg 3 ) < /i 3 , meas(Q g 2 ) < t l 2 with P 3 arbitrary small. @ B 3 and S 3 2 are the sets of near resonant 

3 waves and 2 waves. Here || • || a designates the norm in the Sobolev space H a . 

Theorem 6.1: 

Let B m < Bu < Bm , B m < 1 < Bm ■ Let 0 < zq, < 1 (including iq = v 3 = 0), Bu ©^(^ 2 ), #2 $ ©o 2 . 

Let a > 3/2, a — a > 7, Moa > 0, // 2 ,P 3 < 1- Let. ||U^(0)|| a < Mo a . Let U*(f) be an exact solution of 

the 3D Euler-Boussinesq equations. Let W qg{1) be the solution to the QG equations (3.7), (3.11), (3.12) 
with initial data n^ G U^(0), and w ag(1) the solution to the limit ageostrophic equations (3.17) on the 4 
rays (3.23) with initial data n AG U^(0). Let E(/) be the inertio-gravity waves linear propagator. Then for 
0 < t < Ti 

||U t (i) - W QG (t) - < Ca|(l/,r 2 + al/^/N, for f < N, 

< C Ro a a 3 (l / fi . 2 + 0 ,%/ fi 3 ), for / > N\ (6.3) 

where T\ depends on only on Mo^; C depends only on Mo a , a , Bm (for f < N ) or B m (for / > N). 

Theorem 6.2: 

For n < ^ G U^(f) — W QG {t), under the same conditions as in Theorem 6.1, but with the weaker smoothness 
a — a > 5 we have the estimate, for both inviscid and viscous cases: 

lin^'U^f) - W QG (0||a < CRo a a 3 /p 2 , for / > N 

< Cal/{Np 2 ), for f < N. (6.4) 
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The same estimates hold for ||U (t) — Wq<j(<) — E(— Qw^ G (f)|| a , but with <r — a > 4. Here U designates 
the vertical averaging. 

Remark 6 . 1 : 

For the full error, the above requires smoothness of initial data in Hq. For the convergence of the QG 
component, only H 7. This is a substantial improvement over the H37 in BMN (1996a) and the H is in BMN 
(1996b). The convergence under the Hq smoothness for the vertically averaged fields is rather remarkable, 
as it involves both the QG component and the AG component. It clearly shows that the Dynamical Taylor- 
Proudman theorem established in BMN (1996 a, b) has a modified version for the Bu = 0(1) case, coupling 
QG and AG component. 

In Theorem 6.1, the measures p 2 and p 3 are equal to the measures of the excluded sets of a \ and Bu; 
they reflect estimates of the measure of almost resonant a 2 or Bu. In fact, such estimates imply the “Arnold 
Tongues” of Figure 2: 

Corollary 6.1: 

Let A r > f and 1 /Bu, 0 2 = l/a 2 the resonance parameters, with 1/ Bm < 1 /Bu < 1 / B m , and B m = 0(1). 
Then 

p 3 (l/Bu,9 2 ,0 3 ) = Ca 3, (6.5) 

where C is an absolute constant independent from 9 2 , 9 3 and C = 0(1). The probabilistic (normalized) 
measure of near 3-wave resonances in 1/ Bu on the interval [0, 1 / B m \ is CB m a 3, and the factor a 3 /p 3 within 
the error estimate in Theorem 6.1 reduces to C. 

Corollary 6.2: 

Let N < / and Bu, 0 2 = l/a\ be the resonance parameters, with 0 < Bu < ajj. Then 

fi 3 (Bu, 9 2 ,9 3 ) = 0*0%, (6.6) 

where C* is an absolute constant independent from 9 2 , 9 3 , C* = 0(1). The probabilistic (normalized) 
measure of near 3- wave resonances in Bu on the interval [0,a|] is 0*a|, and the factor a |/// 3 within the 
error estimate in Theorem 6.1 reduces to C* . 

Remark 6 . 2 : 

For N > /, 1 /Bu — ( f/Na 3 ) 2 is the natural resonance parameter; Bu — > +00, N +00 in the strongly 
stratified limit. For N < /; — » 0 as / — > +00 and Bu = (Na 3 / f) 2 is automatically smaller than crj. 

Corollary 6.1 proves that the probability (normalized measure) of near 3-wave resonances at Bu — 
0(1) is very small for the atmosphere on synoptic and mesoscales; H ss 10 km, and even for borderline 
small/mesoscales L k, 50 km, a| « 4 x 10 _i . The 3-wave almost resonant sets for such domains are very 
sparse , quantitatively confirming Bartello’s (1995) picture. 
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Using approximations of solutions of 3 D rotating Euler-Boussinesq equations by solutions of the QG and 
AG equations, we have proven global regularity of solutions of 3 D Euler-Boussinesq equations with arbitrary 
large initial data on arbitrary long time intervals. The global existence theorems in H3 for QG equations 
obtained by Bourgeois and Beale, as well as the conservation of energy on the restricted uncoupled families 
of 4 rays for W ag ensure global existence for W qq + 'Wag- Contrary to Bourgeois and Beale ( 1994 ), we 
do not require any “prepared” initial data: 

Theorem 6.3: 

Let 1 /Bu (/ Gjj 3 ($2 ) 63), O2 (/ €)£ 2 , B m < Bu < Bm ■ Let N > f, v\ — r/ 2 = 0 . Let a > 9 , and M a > 0 , 
T* > 0 be arbitrarily large. Then there exists N* — N*(M a , T* , /j 3 , /j 2 ) such that for |jU^( 0 )|| a < M a and 
N > N * , there exists a unique regular solution U 1 (t) of the 3 D rotating Euler-Boussinesq equations which 
belongs to H a for 0 < t < T* . For M a fixed, T* -4 +00 as N* -* -foo with explicit uniform dependence 
of T* on M a , ^3, /r 2 , N* . Simultaneously, we can take arbitrary large (but bounded) set of initial data: 
M a — > +00 if N* — > +00, for fixed T* . 

If we do not require explicit uniform estimates of T* (equivalently N *) in terms of the geometric param- 
eters <22, <13 and the measures /12, ^3, we can obtain the following “Poor Man’s” long-time existence theorem 
with weaker smoothness assumptions, valid for non-resonant Bu , but holding for all 62 and <7.2 and all 63 
and a . 3 (no constraints on 2-wave resonances ): 

Theorem 6.4: 

Let N > f, u\ — vo — 0 , B m < Bu < Bm, 1 /Bu (/ @3(62,93). Let M4 > 0 , T' > 0 be arbitrary large. 
Then there exists N\ such that if | |U^ ( 0 ) 1 14 < M4 and N > Ni, there exists a. unique regular solution U^(f ), 
0 <t<T of the 3 D rotating Euler-Boussinesq equations, which belongs to H4 as 0 < t < T' . For M4 fixed, 
T' —¥ +00 with Ah — > -foo. Simultaneously we can take arbitrary large (but bounded) sets of initial data: 
M4 +00 if Ni — > +00, for fixed T' . The above holds for all a 2 > «3» 0 < a 2 , a3 < 1 . 

Finally, using techniques from BMN ( 1996 a, b, d) and the fact that the effective spectral eddy viscosities 
lie between v\ and ^ 2 (see Section 3 ), we obtain the regularity for all times for the 3 D rotating Boussinesq- 
Navier Stokes equations (the so called “primitive” equations, not to be confused with equations associated 
to hydrostatic pressure hypothesis). This theorem describes the situation when N is fixed, large enough and 
B m < Bu < Bm- The Navier-Stokes equations are forced by a force F (t) smooth enough. The situation is 
that of non-smooth and arbitrary large initial data in Hq. Then weak Leray solutions exist with maybe a 
blow-up in Hi . 

Theorem 6.5: 

Let vi > 0 , V2 > 0 , N > /, B m < Bu < Bm, 1 /Bu @3(62,63). Let ||F|| 4 + 1 1 F 1 1 1 < Mo 4 for all 
t > 0 , j|Tjt( 0 )||o < Mq. Let T = A-f/ /v m , where u m = min(v i,^ 2 ); T depends only on ||U^( 0 )||o, ^2- 
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Let N' be a number which depends only on M04, v\, v 2, (12 , «3- Then for every fixed N > N' and for any 
weak solution U^(2) of the 3 D rotating Boussinesq-Navier- Stokes equations which is defined on [ 0 ,T] and 
satisfies the classical energy inequality on [ 0 ,T], the following proposition is true: U^(t) can be extended to 
0 < t < 4-00 and it is regular for T <t < +00; it belongs to H4 and [ |U^ (<) | [4 < Ci(Mo4, V\, ^2) for every 
t > T. It F is time-independent there exists a global attractor of the “primitive” equations of geophysics, 
which is bounded in H 4, has a finite fractal dimension, and every weak solution is attracted to the global 
attractor as t — >■ +00. 

Remark 6.3: 

No “preparation” of initial conditions is needed, contrary to the special restricted results of Chemin ( 1995 ). 
This theorem holds for all <22, 0 < 0,2 < 1 . This theorem resolves problems of existence of attractors of for 
the “primitive” equations of geophysics raised by Lions, Temam and Wang ( 1994 ). These authors are able 
to give only conditional theorems assuming a priori regularity in 3 D. We demonstrate the regularizing effects 
of mixed rotation/stratification in the Bu — 0 ( 1 ) regimes. 

7 Concluding remarks 

In this paper, we have treated the problem of strongly stratified limit of rotating Boussinesq equations via two 
complementary approaches. First, we have illustrated the procedure for obtaining the dynamical decoupling 
of the vortical and inertial-gravity wave components of the total flow field. The ‘split’ of the energy transfer 
of the vortical and the wave components is established. As a result, analysis of 3 D Boussinesq system has 
been reduced to several simple, well defined steps. We give explicit phase formulas for the Doppler phase shift 
induced by the interaction between the waves and the mean flow. In the regime of strong stratification and 
weak rotation we proposed regularized reduced equations with control of vertical shearing for all times. The 
regularizing effect of weak rotation principally differs from previously known regularizations which are based 
on vertical viscosity. Second, we have utilized the fact that the time scales of vortical and inertial-gravity 
wave components are disparate. Several useful results can be deduced without going through the steps 
outlined above. For ageost.rophic dynamics, we demonstrate gradual unfreezing of energy cascades as Bu 
varies from zero to infinity. The energy spectrum and the anisotropic spectral eddy viscosity can be derived 
with an explicit dependence on the anisotropic rotation/stratification time scale. This time scale, in turn, 
depends on the aspect ratio parameter (ratio of the vertical and the horizontal length scale). Our analysis 
provides a potential context for investigations of the Garrett-Munk spectrum induced by the interaction 
between inertio-gravity waves and turbulence (Garrett and Munk, 1979 ). 


37 



Acknowledgements 

The authors wish to thank for their support the AFOSR (grant F49620-96-0-0165) and the ASU Centers 
for Environmental Fluid Dynamics and Systems Science and Engineering. The hospitality of the Newton 
Institute (Cambridge) under the special program “Mathematics of Atmosphere and Ocean Dynamics” is 
gratefully acknowledged. We would like to thank Professors C.W. Bardos, P. Bartello, M.J.P. Cullen, M. 
Farge, J.R. Herring, J.C.R. Hunt, J.C. McWilliams, P.S. Marcus, H.K. Moffatt, J. Tribbia, G.K. Vallis, T. 
Warn and V. Zeitlin for very illuminating discussions. 


A Appendix: Dependence of the coefficients Q^n 011 f l = f /N 


The coefficients QjT (/./.) given by Eqs. (2.25)-(2.28) can be expanded in powers of p. =: f/N as follows 

lm£n{k A 77 
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